On the Determination of a Function from an Elliptical Radon 

Transform 



Sunghwan Moon 

Department of Mathematics, Texas A&M University, College Station, Tx77843-3368 



Abstract 

In recent years, Radon type transforms that integrate functions over various sets of el- 
lipses/ellipsoids have been considered in SAR, ultrasound reflection tomography, and radio 
tomography. In this paper, we consider the transform that integrates a given function in 
over a set of solid ellipsoids of rotation with a fixed eccentricity and foci restricted to a hy- 
perplane. Inversion formulas are obtained for appropriate classes of functions that are even 
with respect to the hyperplane. Stability estimates, range description, and local uniqueness 
results are also provided. 
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1. Introduction 

Radon type transforms that integrate functions over various sets of ellipses/ellipsoids 
have been arising in the recent decade, due to studies in synthetic aperture radar (SAR) 
0, 1 i lif, ultrasound reflection tomography 0, [2I, and radio tomography [igI, [H, [13] • In 



particular, radio tomography is a new imaging method, which uses a wireless network of 
radio transmitters and receivers to image the distribution of attenuation within the network. 
The usage of radio frequencies brings in significant non-line-of-sight propagation, since waves 
propagate along many paths from a transmitter to a receiver. Given a transmitter and a 
receiver, wave paths observed for a given duration are all contained in an ellipsoid with foci 



at these two devices. It was thus suggested in [16|, [ITI, ll5| to approximate the obtained signal 



by the volume integral of the attenuation over this ellipsoid, which is the model we study in 
this article. 

Due to these applications, there have been several papers devoted to such "elliptical 
Radon transform." The family of ellipses with one focus fixed at the origin and the other 
one moving along a given line was considered in joj. In the samepaper, the family of ellipses 
with a fixed focal distance was also studied. The authors of |2|, |7| dealt with the case of 
circular acquisition, when the two foci of ellipses with a given focal distance are located on 
a given circle. A family of ellipses with two moving foci was also handled in j^]. 

In all these works, however, the ellipses have varying eccentricity. Also, their data were 

the line integrals of the function over ellipses rather than area integrals. The radio tomogra- 
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phy application makes it interesting to consider integrals over solid ellipsoids. In this article, 
we consider the volume integrals of an unknown attenuation function over the family of el- 
lipsoids of rotation in M" with a fixed eccentricity and two foci located in a given hyperplane. 
We thus reserve the name elliptical Radon transform ReJ for the volume integral of a 
function / over this family of ellipsoids. 

The volume integral of a function /(x) over an ellipsoid of the described type is equal 
to zero if the function is odd with respect to the chosen hyperplane. If the hyperplane 
is given by x„ = 0, we thus assume the function f{x) : M" — M to be even w.r.t. x„: 
f{x',Xn) = f{x', —Xn) where x = (x',x„) G R"~^ x M. 

Given a Radon type transform, one is usually interested, among others, in the following 
questions: uniqueness of reconstruction, inversion formulas and algorithms, stability esti- 



mates, and range description [11|, Il2|. These are the issues we address below. 

The problem is stated precisely in section [21 Various inversion formulas are presented in 
sections [3] and m Sobolev estimates and range description are handled in section O Section [H] 
is devoted uniqueness for a local data problem. 



2. Formulation of the problem 

We consider all solid ellipsoids of rotation in with a fixed eccentricity 1/A, where 
A > 1 and foci located in the hyperplane x„ = 0. We will identify this hyperplane with 
M"^^. The set of such ellipsoids depends upon 2n — 2 parameters, which is n — 2 too many. 
To reduce the overdeterminacy, we require that the focal axis is parallel to a given line, for 
instance, the xi coordinate axis. 

Let u G M""^ be the center of such an ellipsoid and let t > be the half of the focal 
distance. We denote this ellipsoid by Eu^t- Then, the foci are 

ci = (^1 + t, U2, ■ ■ ■ , Un-i, 0) and Ca = (ui -t,U2, - ■■ , u^-i, 0) 

and the points x G Eu^t are described as follows: 

{xi-uiY {X2 - U2)^ xl ^ 2 

A2 A2 - 1 ^ ■ ■ ■ ^ A2 - 1 - ■ 

To shorten the formulas, we are going to use the following notation: 

jj ■= VP^. 

The elliptical Radon transform Re maps a function f{x) into its integrals over the solid 
ellipsoids Eu,t for all u G M""^ and t > 0: 

REf{u,t) = J f{x)dx. 

Eu,t 



Our goals are to reconstruct / from REf and to study properties of this transform. 
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3. Inversion of the elliptical Radon Transform 

In this section, we assume / G C^(]R"). Let xs denote the characteristic function of a 
set S C M'^: 

. . ^ f 1, if X G 5*, 

\ 0, otherwise ' 

Then the elliptical Radon transform can be written as 

REf{u,t) = / XEu.tf{x)dx = C{X) / X\x\<tfi>^Xi+Ui,UX + u,uxn)dx 



t 

n-1 



= C{\) j r'' J }\Xryi + ui.vry + u,ruyn)da{y)dr, 



where u = (ui^u) G M"'"^ x = G M", C(A) = Az/"~^ is the Jacobian factor, and 

a{y) is the surface measure on S"'~^. 

Formula ([1]) can be simplified by differentiation with respect to t and division by 
which yield 



1 d 



REf{u,t) =C{\) j f{\tyi+ui,uty + u,tuyn)da{y) 

dy' 



dt 

bl=i (2) 



C{\) / f{u + {t\y,,tvy),tv^l-\y?) ,^ , 

1 - \u 



\y'\<i 



'12 ■ 



where y' = (yi, y) G M*^"^. 

It is easy to check that -R^; is invariant under the shift with respect to the first n — 1 
variables. That is, if fa{x) := f{x' + a, x„) for x = {x', x„) G M" and a G M"^"*^, we have 

{REfa){u,t) = {REf)iu + a,t). 

Thus, application of the [n — l)-dimensional Fourier Transform with respect to the center u 
seems reasonable. Doing this and changing the variable y' G W^~^ to the polar coordinates 
{e,s) G ^"-1 X [0,oo), we get 

1 



where ReJ and / are Fourier Transforms of ReJ and / with respect to the first n — 1 
coordinates x' of x and u of {u,t), respectively, and 9 = {0i, 9) G S*""^. 
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To compute the inner integral, we use the known identity [sl, ^ 



e^«-^cie = (27r)"/2|e|(2-")/V(„_2)/2(|el)- 



We thus get 



where w„ = 2(27r)("-i)/2C(A). 

We define the back projection operator R*^ for g{u,t) G C^(M"~^ x M+) as 



In fact, R%g{x) is the dual transform not to REf{u,t), but rather to ^REf{u,t), i.e., 

oo 

—REf{u,t)g{u,t)dudt= / f{x)R*Eg{x)dx. 

M"-i R" 

We thus get analogue of the Fourier slice formula. 

Theorem 1. For a function f G C^(M") which is even with respect to x„ the following 
formula holds: 

fiO ^ [^^F^^di^^^ ) 

where Tf is the n- dimensional Fourier transform of f . 

Proof. Let us denote the radial Fourier transform by Hnf{p), i.e., 

oo 

HJXp) :=p'-^/' [ r^l'j^^_^y,{rp)f{r)dr. 



Taking this transform of j^-^^RE f as a function of t, we have for ^ = (^i, ^) = ^' G M" ^, 



oo 1 

3-n 




(5) 

^ \ s-^-^dsdt 



OOnp'-"-'^ I lt-^J^{tp){s\e\)^J^{ts\e\)f{^,^,tuVT^ ^ 



, 
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It is known in js], p. 59 (18) vol.2 or for = 2, p. 55 (35) vol.1] that for a > 0,/3 > 0, 
and, /X > z/ > —1, 

^-M^.+1/2^-M+-+1(^2_^2)1/2m-1/2.-1/2j^_^_^(^(^2_ ^2)1/2) if Q < y < a, 

if a < y < oo. 

To use the above identity, we make the change of variables (s,t) — )■ (x,/3), where t = 
•\/a;2 + /32 and s = x/ y^x^T^ in ([S]), which gives 

oo oo / ~ \ 

^np'^lef- / / |x|J^(p(a;2 + /?2)|)(^2^^2)-^j_^^^|^.|)^~[^^i^^^y^3^j 



dxd(3 







on/2+l„n/2 2-n C C f \ 

C[\) pr^2 - y / p^'/3Hcos(v/7^W/3M/3 if|ri<P, 



otherwise. 

(6) 

Substituting p = |^| yields 

Since / is even in Xn-, the last integral is Fourier transform of / with respect to a; — n, so we 
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Taking Fourier transform of R*^g with respect to x yields 



' 9\u,\ r \^ ^ + ' , ^ ' + ^ I dudx 



A2 



}n-l 



C(A) y e"™ V^'^; y e-^^-^^(M, \x\)dxdu 



(8) 



= (27r)-"/2C(A)i/„^|^|,^,|elj. 

where x = = G M",m = G M"-i and C = (e',^ = (6,e,en) e W\ 

Combining ^ and (E]), we get (g]). □ 

Corollary 2. T/ieorem[I] /ea(is naturally to a Fourier type inversion formula for even func- 
tion, if one supplements (jl]) with the inverse Fourier transform. 

Remark 3. Since = ^^/(O? '"^^ 9^^ for any a,(3 > 0, 

One can also obtain a useful relation with convolution. 

Proposition 4. Let G x [0,oo)), i) = Wcj), f G C^{W) and g = ^^iRsf ■ 

Then we have 

(27r)"/2 d 



where 



oo 



^*0(M,t) = y y g{u-u',t-t')^{u',t')du'dt'. 

M"-i 
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Proof. Note that since if„ is also a radial Fourier transform, we get Hn{f * g) = HnfHnQ 
and HnHnf = f. Taking Fourier transform on g * (p with respect to u and with respect 
to t, we get 

In last line we used equation ([8]). Equation ([7]) implies 

C(A)|en| ^*'^l^A'i.' 

= ^xAr "l^^ V l,A'^'"^"^ 

which proves our assertion. □ 

4. Another Inversion 

In this section, we provide a different inversion formula for the elliptical Radon transform. 
To get it, we start from formula (|2]). Let us define 

00 

G{u,w) ■= j ^^REf{u,t)e'^'"dt 




C(A) J j fiXty^ + Ml, uty + u, tiyyn)e''"'' da{y)dt 



where in the last equality we switched from polar to Cartesian coordinates. 
Theorem 5. Let f G C^(]R") with f{x',Xn) = f{x', — x„). Then we have 



for Xn > 0, where C(A) = Az/"~"^, as before. 



fix) = — / e"*^e'"'^"^'"^e ' ' V^" ' ' G ( ^ — 1 dad-y 

^ ^ ' (27r)"C(A)2 ii,/ V 27 ' 27 ' ^ ' 



Proof. Making the change of variable xi = Xtyi + ui, x = uty + u, x„, = tuyn, we get 



G{u,w) = / f{x)e 



9 I - i' 

e A^^e 



/(x)e 



2 , ,9 9 

TT"l 

1/^ 



where x = and m = G M" ^. Next, make the change of variables 



Xi 



X 



•^l,temp ; •^temp ; 

A U 



and r 



2 I ~ 1 2 2 

I I 7X 

i^^ z/2 ' 



so that 



•^1 •^Ijtemp-^; •^temp^j and X^ U T ■^l^temp \'^temp\ 

The Jacobian of this transformation is 



J 



A 






V 








so that 



2-\ /r — Xij-^j^p — \Xtemp\'^ 



dx 



C(A) 



-dxidxdr. 



\temp l^temp 

Let the function k{x,x,r) be defined by 

/(Axi, z/x, z^A/r — x\ — |xp) 







2^r — x\ — |x'^ 



< |xip + |xp < r, 
otherwise. 



Since / is even in x„, it is sufficient to consider the positive root of y P. Then we 

can rewrite G{u, w) as 



4 ■ ifip 

— ^IW- — - 



G{u,w) =C(A)e™i*e™^ / k{xi,x,r)e''"'e-''''^e-^''^ dxidxdr 



4 ■ isp 



C(A)e™i*e™1^ A' ( 2^, 2—, -w 

A 



where for a = {ai,a') G M x 



bn—2 



K{a,-f) = I e-'''<'''^^^e-''"'k{xi,x,r)dxidxdr 

i jHi^ f —ai\ ~a'u \ 
e G{ ——, -— , -7 . 



C(A) V 27 27 

Since k{x, x, r) is 

—V f e"''=''e"''-^e'^'K{a,-f)dad-f 

(27r)'" jRn 

we get for > 0, 

(~ 2 I ~ 1 2 2 



/ ■.Vir(a,7)dad7 (9) 



(27r)«C(A) 7r 

(27r)«C(A)2 A„ ' ^^27'27'^'''"''^- 



□ 



5. A Sobolev Estimate and Range description 

In this section we obtain a Sobolev estimate for the elhptical Radon transform and 
describe its range. Let consist of all L^-functions even in the x„ variable. If is a 

function on MJ^~^ x [0, 00), we define by 







V^^'dtdu. 



Note that || ■ ||i„„ is norm on C^{W~^ x [0,oo)). Also, let Hl_^{W^^ x [0,oo)) = {g : 
< 00}/ ~, where the equivalence relation ~ /i is defined by g{u,t) = h{u,t) + a{u) 
for some function a{u). Then ifj^_„(R"^^ x [0, oc)) is a vector space and || ■ ||i_„ becomes a 
norm on this space. Let EC^(K"-~^ x [0, 00)) be defined as the quotient space C°°(M"'~^ x 
[0,00))/ ~ and [g] G EC^{W~^ x [0,oo)) for g G C~(M"-i x [0,oo)) be the equivalence 
class of g. 

Theorem 6. Hl_^(R"'~^ x [0, 00)) is complete with respect to the norm \ \ ■ \ \i-n (ind has 
EC^iW"-^ X [0,00)) 05 a dense subset. 
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Proof. To show this space is complete with respect to || ■ ||i-n, let C Hl_^(W' 

[0,00)) be a Cauchy sequence. Then {■§igi{u,t)}'^i is a Cauchy sequence in L 



n-l 



2mn-l 



[0, 00), t with a measure t ^dtdu. Hence, there exists h in L 



X [0, 00), t 



l-n\ 



such 



that 



Define a function (7 by 
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—gi{u,t) - h{u,t) 



f^^'dtdu 0. 



g{u, t) := / r)dr. 



Clearly, -§igiu,t) is /i(u,t), [g] is in iJ|_„(M""i x [0,oo)), and \\[gi] - [g]\\i^n goes to zero. 
This proves the completeness of Hl_j^{M."'~^ x [0, 00)). 

To show that EC^{W-^ x [0, 00)) is dense in Hl_^{W'-^ x [0, 00)) with respect to 1 1 ■ 1 
pick [g] e i/^„„(M"-i X [0,oo)). Then §-^g{u,t) is in L2(M"-i x [0, cx)), t^""), so there exists 
a sequence in C^{W-^^ x [0, cx))) such that 



d 

hi{u,t) - Q^g{u,t) 



t^-^'dtdu 0. 







For each i, let gi{u,t) = hi{u,r)dr. Then the sequence {gi}'fZi ^ C*i 



00 /"Iipn— 1 



X [0,00)) 



satisfies hi{u,t) = ■^gi{u,t) and 11(71 — 5'||i-n goes to zero. Hence, for the sequence {[s'j]}^! E 



EC: 



00 /"Iipn—l 



X [0,00)), \\[gi\ - [g]\\i-n goes to zero. 



□ 



We get a following estimate and range description. From now on, let the equivalent class 
[g] be denoted by g. 



Theorem 7. For f E L 



2^'^^)n^ 



is equal to | li-n? up to a constant multiple factor. 



Furthermore, Re is an isomorphism between and Hl_^{M."' ^ x [0,oo)). 

Proof. Let g = Rsf ■ The Plancherel formula implies 



1 d 



dt 



du 



1 d 

Hn^,-^^g{v.-){p) 



drj 



1 d 



(P) 



(10) 



dr]. 



Note that f{x) = HnHnf{x) for a radial function / on M". By the Plancherel formula, we 

get 
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oo oo 

I |/(r)|V-Mr = I \HJ{p)\'p-~Hp. 





Combining two above equations f fTOj) and ffTTj) . we get 



'11^ 



li/ 1 ll-n 



-1 



1 9 



dt 



-1 



1 d 



R"-i 
2n/2+l^n/2 



Ji-i 



1 d 
t^dt 

P 
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A u 



ip) 



p" dpdr] 



\V\ 



^ ' A'z/' u 



R"-i \ri\ 

2"/2+i7r"/2c(A) / \f{r],p)\^dridp. 



dpdf] 



In the fourth hues, we used equations (|6]) and ((Tj), and in the last hne, we changed to a variable 
Pnew = Vp2 - Ivl^/i" and then variables r]i^new = Vi/^, Vi,new = r]i/u for i = 2, ■■ ■ ,n-l. 
To show that Re is an isomorphism, pick g G £'C^(M"'"^ x [0, oo)). It is routine to check 
e C~(M'^). / e L2(R'^) n ^^(M") follows from equation ©, so / G ^^(M"). 
Therefore, .R^; is an isomorphism since C^(]R") fl Lg(M'^) and EC^{M.'^~^ x [0, oo)) are dense 
in Lg(M") and Hl_.^{W^^ x [0, oo)) with respect to || ■ II2 and || ■ ||i-n, respectively. □ 



6. Uniqueness for the Local Problem 

Theorem [7] implies / G Lg(M") is uniquely determined by Re/- The question arises if / 
is uniquely determined by some partial information. The approach in this section is similar 
to the one in 

Theorem 8. Let vP E W^~^, e > 0, and R > be arbitrary. Let f G Lg(]R") and suppose 
g = ReJ is equal to zero on the open set 

UT,e = {{u, t) G M""^ X [0, 00) : |n - n°| < e, < t < T}. 

Then f equals zero on the open set 



Here x = {xi,x,Xn) G and n° = {ui,u^) G M" ^. yl/so, is equal to zero on the open 
cone 

Wt = {{u, t) e M"-^ X [0, T) : |n - uol + 1 < T}. 

Proof. Without loss of generality, we may assume = 0. Let / G C°°(M"). Clearly, g is 
also differentiable. Differentiating REf{u,t) with respect to yields 

d f d 

—REf{u,t) = C(A) / X\x\<t^f{>^Xi + Ui^VX + u,vxn)dx 
= C(A)- / — /(Axi + Ui, + £t, 

\x\=t 

Here we used equation ([1]) and the divergence theorem. Using equation ([2]), we get 
d f 

—RE{xif){u,t) =C{\) / {uXi + Ui)f{Xxi + ui,ux + u,uXn)da{x) 

\x\=t 

f d d \ 

= C{\) \ tiy—g{u, t) + Ui^g{u, t) 1 . 

Let the linear operator Di be defined by the above expression. Then -^RE{xif){u,t) is 
Dig(u,t). By iteration, we obtain -^RE{p{x')f) = p{D)g where p is an n — 1-variable 
polynomial. If g is zero in UT,e, then p{D)g is also zero in UT,e- Then we have for any point 

{U,t) G UT,e, 

d f 

:RE{p{x')f){u,t) = C(A) / p{\xi + Ui,vx + u)f{\xi + Ui,iJX + u,vXn)da{x) 



dt 

\x\=t 

= C{\) [ p{u + (Ayi, uy))f{u + {Xy,, uy), v^t'' - \y\^)—^== = 0. 

J - \y\ 

\y\<t 

For fixed u and t, choose a sequence of polynomials such that Pi{u + {Xyi, uy)) converge to 
f{u + (A?/i, i/y, v^t'^ — \y^)) uniformly for \y\ < t and y = {yi,y) G MJ^^^. It follows that 
/ = in Vr and that = in Wt- 

To extend to the case / G L^(M"), we choose a sequence of gi G C°°(]R") such that 
\\9i ~ fl'lli-n — ^ where g = ReI- We may assume (^j = in Ux^^ei C f/^.e such that Tj — )■ T 
and ej — i- e. By theorem [3, there exist /j such that REfi = di and fi converge to / in the 
L^(]R") sense. Thus fi is zero on Vr^, so / is also zero on Vr- Then g also vanishes in Wt- □ 

Corollary 9. If g = in a strip {{u,t) G M"~^ x [0, cx)) : |m — Mo| < e, < t < oo} then f 
and g are zero everywhere. 



12 



7. Conclusion 



We describe two different ways of determining a function / from its n- dimensional ellipti 



cal Radon transform Re/ arising in radio tomography imaging [16|, [151, Il7[. We also present 
a Sobolev stability estimate, range description and a local uniqueness result for Re- 
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